arXiv:1505.07864vl [astro-ph.GA] 28 May 2015 


Draft version June 1, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


THREE CANDIDATE CLUSTERS AROUND HIGH REDSHIFT RADIO-LOUD SOURCES: MG1 J04426+0202, 

3C 068.2, MS 1426.9+1052 

J.R. Franck 1 , S.S. McGaugh 1 , J.M. Schombert 2 
Draft version June 1, 2015 

ABSTRACT 

We present near-infrared observations of the environments around three radio-loud sources (MG1 
J0442+0202, 3C 068.2, and MS 1426.9+1052) at redshifts z = 1.10,1.57, and 1.83 (respectively), that 
are surrounded by near-infrared galaxy overdensities. Overdensities with respect to field counts were 
found to be significant up to 19a, with twelve times the expected number of galaxies within the inner 
regions of the densest proto-cluster. Color-magnitude relations are constructed in K Sl J — K s , with 
each candidate cluster exhibiting a feature consistent with the beginnings of a red sequence. Galaxy 
models based on the redshift of the radio source are used to compare expected color-magnitude 
relations for a given formation epoch with the observed red sequence of each candidate, and are found 
to be consistent with an old (zf > 5) formation epoch for a few bright, red galaxies on the red 
sequence. 

Subject headings: galaxies: clusters: general - galaxies: clusters: individual(CL J0442+0202, 3C 068.2, 

MS 1426.9+1052) 


1. INTRODUCTION 

The largest overdensities in the early universe are ob¬ 
served today as galaxy clusters. These initial density 
perturbations can be used as probes over a range of 
epochs for the cosmological parameters that character¬ 
ize our universe (Cen & Ostriker 1994; Vikhlinin et al. 
2009). Cluster observables (number density, mass) as a 
function of redshift provide constraints for and U m . 
The mass and formation timescale of high redshift clus¬ 
ters compared with simulations of large-scale structure is 
a powerful test of ACDM (Mortonson et al. 2011). Re¬ 
cent discoveries are challenging the boundaries at how 
rapidly large structures can form {z ~ 2 from Gobat 
et al. 2013; z ~ 3.8 from Lee et al. 2014a; 2 ~ 4 from 
Venemans et al. 2002), and how massive they can become 
(e.g. ‘El Gordo’ M > 10 15 Mq at 2 ~ 0.9 in Menanteau 
et al. 2012). Identification of the epoch of cluster for¬ 
mation is an important test of cosmology and the dark 
energy equation of state. 

Previously these high redshift structures were identi¬ 
fied via galaxy overdensities in near-infrared (NIR) pho¬ 
tometric surveys (Stern et al. 2003; Papovich 2008; Eisen- 
hardt et al. 2008; Papovich et al. 2010). Clusters are also 
associated around radio-loud sources (Hall et al. 2001) 
out to redshifts z ~ 4 (Venemans et al. 2002). Utilizing 
the intra-cluster medium (ICM) as an identifier, diffuse 
X-ray emission has been useful in increasing the num¬ 
ber of known high redshift clusters (Gobat et al. 2011; 
Willis et al. 2013), while the Sunyaev-Zel’Dovich (SZ) ef¬ 
fect is becoming increasingly powerful with the success of 
the South Pole Telescope (Keisler et al. 2011), Atacama 
Cosmology Telescope (Sievers et al. 2013), the Planck 
satellite (Planck Collaboration et al. 2013a) , and with 
the Combined Array for Research in Millimeter-wave As¬ 
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tronomy (CARMA 3 ; Muchovej et al. 2007; Culverhouse 
et al. 2010; Mantz et al. 2014). 

We present three cluster candidates at z > 1, two of 
which have not been identified previously. High-redshift 
radio-loud sources were targeted for NIR photometry to 
search for a galaxy overdensity in the surrounding re¬ 
gion, which has proven to be a robust identifier of struc¬ 
ture (Stern et al. 2003; Haas et al. 2009; Galametz et al. 
2010). Section 2 describes the NIR observations and pho¬ 
tometry, while the overdensity associated with each clus¬ 
ter candidate is discussed in Section 3. The photometry 
of each candidate cluster is presented in Section 4 and 
analyzed in Section 5, with a summary presented in Sec¬ 
tion 6. 

The work presented here assumes a ACDM concor¬ 
dance cosmology, with Ha = 0.7, a matter density of 
fl m = 0.3, and H 0 = 70 km s _1 MpcW 

2. OBSERVATIONS 

Three candidate clusters were observed with Kitt Peak 
National Observatory’s Mayall 4-m telescope equipped 
with the SQIID Camera in the NIR J, H , K s and L pass- 
bands between January 2006 to April 2007. A summary 
of observations can be found in Table 1. Each passband 
is projected onto a 440 x 460 pixel InSb detector (cover¬ 
ing 172" x 179" in RA/DEC). Individual frames were co¬ 
added and corrected for cosmic rays with IRAF 4 (Tody 
1986, 1993) IMCOMBINE routined 

The candidate cluster around 3C 068.2 was also imaged 
with the NEWFIRM camera on the 4-nr in October 2008. 
NEWFIRM is a wide-held camera that spans 28' x 28', 
and images were taken in J, H, and K s Liters. The clus¬ 
ter fell entirely within one quadrant of the NEWFIRM 
camera, so we devised the “four-shooter” observational 

3 http://www.mmarray.org 
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Table 1 

Observations Log 


Date 

Field 

Instrument 

Bands 

Typical Seeing 

7-9 January 2006 

CL J0442+0202 

SQIID 

J,H,K S ,L 

0.5”-0.8” 

7-9 January 2006 

3C 068.2 

SQIID 

J,H , K S ,L 

0.5”-0.8” 

7-9 April 2007 

MS 1426.9+1052 

SQIID 

J,H , K S ,L 

1.1” 

12-13 October 2008 

3C 068.2 

NEWFIRM 

J, H, K s 

1.1”-1.3” 


Note. — Observations Summary for the three selected candidate clusters. 


mode in which the telescope is shifted from quadrant 
to quadrant to maximize time on target while simulta¬ 
neously observing sky. These frames were processed by 
the NEWFIRM Quick-Reduce Pipeline. The galaxy pho¬ 
tometry between the NEWFIRM and SQIID images were 
within the photometric errors, increasing our confidence 
in the calibration and detection procedures. 

Sky subtraction was accomplished by chopping on and 
off field while dithering between frames. The resulting 
sky flats displayed variations of only 0.02% based on 
sky box determination of the frame sky value. Hystere¬ 
sis from a bright star in the I\ s frames for cluster MS 
1426.9+1052 produced some large scale features (visible 
under high contrast, see Figure 5). However, these vari¬ 
ations are less than 0.04% the background and do not 
effect the galaxy photometry routines which used local 
sky. Both J and K s frames were used for galaxy identi¬ 
fication and each object was confirmed visually in both 
frames before inclusion in our analysis. 

Calibration of the stacked images utilized the 2MASS 
Point Source Catalog (PSC) for stars within the individ¬ 
ual frames in the J and K s filters. With the wide field 
of view (FOV) of NEWFIRM , a large sample of stel¬ 
lar sources were used for calibration (474 matches for J 
and 577 stars for K s ), providing photometric solutions 
with intrinsic scatter of aj = 0.05 mag and gks = 0.09 
mag. A photometric solution for the SQIID instrument 
was found in a similar way, but the smaller FOV of each 
image provided fewer 2MASS calibration sources. CL 
J0442+0202 and 3C 068.2 were imaged on the same run 
in January 2006, while MS 1426.9+1052 was observed 
in April 2007. Unfortunately, no usable 2MASS stars 
were in the latter field, so the photometric solutions for 
MS 1426.9+1052 were derived using 2MASS stars within 
other science frames taken during the April 2007 run. 

Source Extractor (SExtractor; Bertin & Arnouts 
1996) was used to create an initial list of sources from 
both J and K s frames. Galaxies and stars were distin¬ 
guished using the stellar index parameter. A few of the 
faintest objects were below the detection threshold in one 
filter (typically J), and were confirmed visually, in both 
filters. A few objects were visible in the I\ s frames, but 
rejected as they were below the threshold or not visible 
in the J frame preventing a color determination. The 
final set of objects are marked in Figures 1, 3 and 5, 
all clearly visible under high contrast. Overdensity esti¬ 
mates are made using a conservative limiting magnitude 
in the K s frame, well below the detection threshold. 

After the photometric solution was applied to the 
candidate cluster galaxies, a color-magnitude relation 
(CMR) was constructed using 4" diameter apertures. 
The aperture size was determined by building a curve 


of growth for the samples and examining the source light 
profiles with IMEXAM. The 4" diameter apertures were 
found to be a suitable balance of light collected versus 
noise levels. 

3. CLUSTER-CANDIDATE GALAXY OVERDENSITIES 

The most basic technique of identifying clusters is by 
finding overdensities of galaxies with respect to field 
counts. The lack of a large overdensity, conversely, does 
not imply that no cluster exists, as cluster richness and 
the system’s mass were not found to be strongly cor¬ 
related (Planck Collaboration et al. 2013b), although 
other recent results have suggested the opposite (Rozo 
& Rykoff 2014). Following the example of Hall & 
Green (1998) and Hall et al. (2001), counts of galaxies 
brighter than K s = 20.5 surrounding radio-loud quasars 
consistently show overdensities relative to field counts 
(13.7 ± 1.5 arcmin~ 2 ). Listed in Table 2 is each radio 
source, its position and redshift, and number of galaxies 
identified within ~ 1' of the radio source. 

For each cluster candidate, a peak overdensity was de¬ 
termined by counting galaxies with K s < 20.5 around 
the radio source, and comparing this to the expected 
field counts. The difference between the counted and 
expected (in units of er = 1.5 counts arcmin -2 ) are pre¬ 
sented in Table 2 as the galaxy overdensity. The peak 
surface density for each candidate is computed as the 
number of galaxies divided by the same area as the peak 
overdensity, written in units of galaxies per square ar- 
cminute. 

The lowest redshift source in this sample was found 
to have the largest overdensity of 19.2er above the mean 
value, followed by 3C 068.2 with a 14.4<r overdensity. The 
greatest redshift candidate cluster (at £ = 1.83) around 
MS 1426.9+1052 had the smallest overdensity of 12.6o\ 
This can be attributed to the off-center position of MS 
1426.9+1052 in the SQIID image (see Fig 5) effectively 
limiting the potential number of galaxies, coupled with 
the smaller amount of time the cluster candidate has had 
to dynamically assemble (~2 Gyrs between 0 = 1.11 —► 
1.83). 

The expected number of field galaxies that reside 
within the region associated with each overdensity can 
be estimated utilizing the counts of Hall & Green (1998). 
Within the areas associated with the peak overdensi¬ 
ties (e.g. > 10cr), the estimated number of galaxies un¬ 
associated with the cluster is ~ 4/13, ~ 4/10, and ~ 4/10 
for CL J0442+0202, 3C 068.2, and MS 1426.9+1052, re¬ 
spectively. As each cluster candidate’s overdensity is con¬ 
tained within approximately the same area on the sky, 
they have approximately the same number of expected 
field galaxies. 

It is important to note that galaxies which exist out- 
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Table 2 

Summary 


Field Redshift RA DEC Galaxies Measured Overdensity Peak Surface Density 

arcmin ~ 2 

CL J0442+0202 z = 1.10 4:42:23.7 +02:02:20 48 19.2<r 42.4 

3C 068.2 z = 1.57 2:34:23.8 +31:34:17 31 14.4cr 35.4 

MS 1426.9+1052 z = 1.83 14:29:24.1 +10:39:15 33 12.5<r 32.4 

Note. — Summary Table for the three selected candidate clusters. All radio-loud sources have significant galaxy 
overdensities surrounding them with respect to the NIR field counts provided by Hall & Green (1998). The Surface 
Density (in galaxies arcrnin -2 ) is measured within the same area as the overdensity for sources K s < 20.5. 


side of the largest overdensity can still be associated with 
the cluster. For instance, out of the 6 spectroscopically 
confirmed members within CL J0442+0202, Galaxy 4 in 
Tab 3 is in an area outside an even modest overdensity 
(< 3cr). Within the cluster discovered by Papovich et al. 
(2010) at z = 1.62, 8 of the 11 spectroscopically con¬ 
firmed members would exist outside of a ~ lex overdense 
region if these galaxies were inserted into our data. This 
suggests that at these redshifts, clusters are still in the 
process of assembly with galaxies infalling. 

By computing the expected number of galaxies for a 
projected area in square arcseconds around the radio 
source, an overdensity factor can also be computed. The 
inner regions show a factor of ~ 11.6, 6.9, and 6.1 times 
the expected number of galaxies near CL J0442+0202, 
3C 068.2, and MS 1426.9+1052, respectively. These over¬ 
densities indicate that the three regions are candidates 
for clusters in the process of assembly. 

4. CLUSTER CANDIDATES 
4.1. CL J0442+0202 

Stern et al. (2003) discovered a NIR galaxy overden¬ 
sity surrounding the radio source MG1 J04426+0202 at a 
redshift of z = 1.11, and did follow-up observations spec¬ 
troscopically and in the X-ray with Chandra . Five other 
galaxies were spectroscopically confirmed to be cluster 
members, and a number of X-ray point sources were re¬ 
solved. However, Stern et al. (2003) were not able to 
detect diffuse X-ray emission from the ICM, only from 
the central radio source and other point sources. The 
NIR photometry of the cluster was measured in the I- 
band and Lick K' to a limiting magnitude of K' < 20. 

We present SQIID photometry to a deeper limiting 
magnitude of K s < 20.5. A comparison of 10 galax¬ 
ies with K' < 19.1 magnitudes in the sample of Stern 
et al. (2003) with ours reveals a systematic offset of 
K' — K s = 0.32+0.16 mag. This provides a basic check of 
the photometric solution for the SQIID camera, as well as 
a comparison to published work. The offset could be at 
least partially attributed to the blueward shift of the K' 
filter. The uncertainty in the offset approximately mir¬ 
rors the uncertainty within our own photometry, while 
the uncertainty of the K' photometry from Stern et al. 
(2003) is unknown. 

Within the field of view (shown in Fig 1), we present 
the photometry of 48 galaxies in Table 3 and the resulting 
color-magnitude relation in Fig 2. Spectroscopically con¬ 
firmed members (Stern et al. 2003) are consistent with 
being on a red sequence, which may indicate the initial 
establishment of the prominent red sequences of early- 
type galaxies common in rich, low-redsliift clusters (e.g. 


Table 3 

CL J0442 Sources 


Galaxy 

ARA 

(”) 

ADEC 

(”) 

Ar 

(”) 

K s 

(mag) 

J — K s 

(mag) 

i 

0.000 

0.000 

0.000 

17.06 

± 

0.19 

1.73 

± 

0.21 

2 

27.571 

-37.814 

46.798 

18.11 

± 

0.20 

1.52 

± 

0.22 

3 b 

42.527 

-17.161 

45.859 

17.95 

± 

0.20 

2.15 

± 

0.23 

4 a 

28.336 

-15.598 

32.345 

18.67 

± 

0.20 

1.73 

± 

0.24 

5 

4.610 

-13.686 

14.441 

20.70 

± 

0.49 

0.79 

± 

0.58 

6 C 

-9.483 

-11.521 

14.922 

18.12 

± 

0.20 

1.30 

± 

0.22 

7 b 

-23.240 

-9.315 

25.037 

17.59 

± 

0.19 

1.73 

± 

0.21 

8 

6.097 

-7.555 

9.709 

17.75 

± 

0.19 

0.88 

± 

0.21 

9 

23.292 

-3.630 

23.574 

18.80 

± 

0.21 

1.25 

± 

0.23 

10 

-0.907 

4.254 

4.350 

19.54 

± 

0.26 

1.74 

± 

0.37 

ll a 

-7.788 

5.418 

9.487 

18.11 

± 

0.20 

1.63 

± 

0.22 

12 c 

-2.979 

11.825 

12.194 

19.17 

± 

0.22 

0.86 

± 

0.25 

13 

13.204 

14.311 

19.472 

18.50 

± 

0.20 

1.63 

± 

0.23 

14 a 

3.693 

14.121 

14.596 

19.20 

± 

0.22 

1.59 

± 

0.28 

15 

-28.050 

29.989 

41.063 

19.07 

± 

0.21 

1.59 

± 

0.26 

16 

-21.010 

-46.483 

51.011 

19.51 

± 

0.25 

0.72 

± 

0.28 

17 

15.502 

34.860 

38.151 

19.30 

± 

0.23 

1.97 

± 

0.33 

18 b 

17.970 

40.284 

44.111 

18.38 

± 

0.20 

1.62 

± 

0.23 

19 

33.668 

28.809 

44.311 

19.74 

± 

0.26 

1.93 

± 

0.40 

20 

-12.425 

-8.076 

14.819 

19.53 

± 

0.23 

1.95 

± 

0.38 

21 

-8.792 

0.011 

8.792 

20.16 

± 

0.41 

1.39 

± 

0.58 

22 

-31.158 

-2.055 

31.226 

20.04 

± 

0.30 

0.73 

± 

0.34 

23 

-33.491 

1.069 

33.508 

19.62 

± 

0.25 

1.17 

± 

0.30 

24 

-11.603 

-34.214 

36.128 

19.25 

± 

0.22 

1.86 

± 

0.31 

25 

-32.993 

-37.985 

50.313 

18.05 

± 

0.20 

0.72 

± 

0.21 

26 

-15.997 

-49.336 

51.865 

19.34 

± 

0.23 

1.38 

± 

0.28 

27 

66.671 

-35.440 

75.505 

18.99 

± 

0.21 

1.08 

± 

0.24 

28 a 

4.659 

23.257 

23.719 

19.12 

± 

0.22 

1.82 

± 

0.29 

29 

-6.299 

26.279 

27.024 

19.10 

± 

0.21 

1.89 

± 

0.28 

30 

20.863 

29.580 

36.197 

19.92 

± 

0.30 

0.91 

± 

0.34 

31 

36.724 

19.865 

41.752 

20.09 

± 

0.31 

1.47 

± 

0.43 

32 

11.888 

-43.960 

45.539 

19.24 

± 

0.22 

1.72 

± 

0.31 

33 

56.662 

20.766 

60.347 

19.67 

± 

0.27 

0.99 

± 

0.31 

34 

40.150 

49.944 

64.082 

20.14 

± 

0.34 

0.85 

± 

0.40 

35 

43.020 

55.760 

70.426 

19.22 

± 

0.23 

1.13 

± 

0.26 

36 

23.452 

56.852 

61.500 

19.18 

± 

0.22 

1.45 

± 

0.27 

37 

51.838 

50.321 

72.245 

18.54 

± 

0.20 

2.58 

± 

0.30 

38 

25.158 

-0.324 

25.160 

19.53 

± 

0.24 

1.59 

± 

0.32 

39 c 

7.734 

52.856 

53.419 

18.55 

± 

0.20 

1.71 

± 

0.24 

40 

14.390 

38.199 

40.820 

20.47 

± 

0.38 

0.87 

± 

0.44 

41 

68.910 

39.558 

79.457 

19.06 

± 

0.22 

2.02 

± 

0.33 

42 

42.107 

48.183 

63.989 

19.64 

± 

0.25 

1.49 

± 

0.33 

43 

23.571 

26.444 

35.424 

20.29 

± 

0.34 

0.84 

± 

0.40 

44 

32.678 

10.572 

34.346 

19.87 

± 

0.27 

2.35 

± 

0.60 

45 

-29.289 

13.467 

32.236 

19.49 

± 

0.24 

2.20 

± 

0.41 

46 

-27.731 

-15.182 

31.615 

20.02 

± 

0.30 

1.90 

± 

0.52 

47 

25.050 

-16.000 

29.724 

19.86 

± 

0.27 

1.31 

± 

0.35 

48 a 

-2.0272 

-2.370 

3.118 

18.43 

± 

0.20 

1.76 

± 

0.24 


Note. — The central radio source MG1 J04426+0202 (Galaxy 1) 
and the sources around it within the SQIID field. Galaxy 1 is located 
at RA=4:42:23.7 DEC=+02:02:20 with separations (ARA,ADEC) 
computed with respect to it. 

a Spectroscopically confirmed cluster members (Stern et al. 2003). 
b Non-members in the field (Stern et al. 2003). 

c X-ray point sources found with Chandra that are not cluster mem¬ 
bers (Stern et al. 2003). 
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Figure 1. The SQIID K s image of CL J0442+0202 that spans 120” X 120”. The radio source MG1 J04426+0202 is indicated by a red 
arrow and square, while spectroscopically confirmed cluster members (Stern et al. 2003) are inside magenta ellipses. The green boxes are 
the remainder of the sources measured photometrically. 


Coma; Eisenhardt et al. 2007). There is also a bright 
interloper (K s = 17.59) posing as a red sequence mem¬ 
ber. There are also foreground and background AGN 
(bright X-ray point sources found with Chandra ) that 
are entangled within the CMR. The overdensity remains 
significant, but the aforementioned interlopers urge cau¬ 
tion in relying on purely photometric results. 

4.2. 3C 068.2 

The photometry for the galaxies around 3C 068.2 are in 
Table 4. Within ~1 arcminute of the central radio source 
(Galaxy 1 in the table), which corresponds to roughly 
~ 0.5 Mpc in the assumed cosmology (Wright 2006), 31 
extended sources were identified. The selected sources 
are shown in Fig 3. 

A K 3 , J—K s color-magnitude relation of the sources is 
in Fig 4, displaying a flat red sequence (RS) of J—K s « 2. 
The two bright galaxies on the RS with K s = 16.30,17.55 
likely contain AGN, as they were identified as X-ray point 
sources by Wilkes et al. (2013) with Chandra . 

It is interesting to note that the central 3C source in 
this case is neither on the red sequence, nor is it the 
brightest cluster galaxy. Instead, it is found in the so- 
called blue cloud, inhabited primarily by star-forming 
galaxies. Barthel et al. (2012) surveyed a number of 
3C sources with Herschel and integrated the far-infrared 
(FIR) dust emission to calculate a star formation rate 
(SFR) of each galaxy. 3C 068.2 has a calculated SFR of 
~ 390 Mq yr -1 . This exceptionally high SFR cannot be 
sustained for a Hubble time without exceeding the total 
mass of the largest known galaxies. This might suggest 


that this galaxy has been caught in the act of forma¬ 
tion prior to its transition from the blue cloud to the red 
sequence. 

4.3. MS 1426.9+1052 

MS 1426 is the most distant candidate cluster in this 
sample, and it contains a large number of potential clus¬ 
ter galaxies, with 34 measured photometrically in Ta¬ 
ble 5. The SQIID image (Fig 5) surrounding the radio 
source was slightly off-center, so the selected galaxies are 
not evenly distributed from the radio source. The color- 
magnitude relation is shown in Fig 6, with a RS color 
consistent with the CMRs for CL J0442+0202 (Fig 2) 
and surrounding 3C 068.2 (Fig 4). 

Richards et al. (2009) computed photometric red- 
shifts of regions surrounding radio-loud sources within 
the SDSS database. There are two galaxies near MS 
1426.9+1052 with photometric redshifts of z = 1.83, but 
were unfortunately off of the edge of the SQIID image. 
However, this provides a further indication that the field 
hosts a possible structure, in addition to the presence of 
a red sequence in the color-magnitude relation and the 
associated galaxy overdensity that was discussed in Sec¬ 
tion 3. 

5. DISCUSSION 

Perhaps one of the oldest galaxy correlations is the 
one between galaxy absolute luminosity and global color, 
the color-magnitude relation (i.e., the CMR; Baum 1959; 
Faber 1973; Visvanathan & Sandage 1977). The CMR 
has been demonstrated to exist from the near-UV (Kavi- 
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Figure 2. The radio-loud source (MG1 J04426+0202) at the center of CL J0442+0202 is shown by a red square as the brightest cluster 
galaxy (BCG), and is also a part of the red sequence (J — K s = 1.73) in this color-magnitude diagram. The next brightest galaxy along the 
red sequence (a green triangle at K s = 17.59, J — K s = 1.73), is a spectroscopically confirmed interloper. Magenta squares are confirmed 
members (Stern et al. 2003), while foreground/background X-ray emitting AGN are represented by blue downward triangles. Galaxies 
that are found within the region corresponding to at least a 10(7 galaxy overdensity have the largest symbols, while the smallest points 
correspond to galaxies located in a region with an overdensity < 3a. Within the largest overdensity region, approximately 4 sources are 
expected to be unassociated with the cluster. 


raj et al. 2007) to near-IR wavelengths (Chang et al. 
2006). A majority of these CMR studies can be sum¬ 
marized such that (1) the correlation is one of redder 
galaxy colors with increasing galaxy luminosity, (2) the 
correlation exists over a range of galaxy environments, 
but is strongest in the densest galaxy regions (Cooper 
et al. 2008), and (3) there exist distinct red and blue 
components to the CMR, with the red component being 
composed of early-type galaxies (ellipticals and SO’s, see 
Bower et al. 1992; van Dokkum et al. 1998; Hogg et al. 
2004; Lopez-Cruz et al. 2004; Baldry et al. 2004; Bell 
et al. 2004; Bernardi et al. 2005; McIntosh et al. 2005; 
De Lucia et al. 2006; Gallazzi et al. 2006). 

The CMR is a member of a number of fundamental 
relations that relate the global properties in early-type 
galaxies, such as mass and kinematics (Dressier et al. 
1987). Various formation and evolutionary scenarios 
are proposed to explain these fundamental relationships 
(Gavazzi et al. 1996). For example, galaxy interactions 
can remove the gas used for star formation and produce 
a relationship between age and galaxy mass (Moore et al. 
1999). Ram pressure stripping can also strip galactic gas 
and halt star formation (Tecce et al. 2010). However, 
the comparison between spectroscopic indices and galaxy 
color indicates that the CMR is most strongly linked to 
the chemical enrichment history of a galaxy, as redder 
colors map directly into increasing metallicities (Thomas 
et al. 2005; Bernardi et al. 2005). The slope and scatter 


of the CMR can be reproduced if the redshift of galaxy 
formation is greater than 5 (Bower et al. 1998) even with 
a small amount of mass growth by mergers or bluer colors 
with recent star formation. 

With respect to cluster populations, the CMR has also 
been used to define the ‘red envelope’ or red sequence (El¬ 
lis 1988), the boundary between old, quiescent galaxies 
and star-forming systems. The small amount of scatter 
for the red sequence of the CMR has been used to ar¬ 
gue that early-type galaxies are coeval and have evolved 
passively since their formation epoch (McIntosh et al. 
2005; Schombert & Rakos 2009b). The blue population 
(the ‘blue cloud’) can be compared to the red population 
by assigning the fraction of galaxies below 0.2 magni¬ 
tudes of the red sequence; this fraction (/s) becomes 
an indicator of the global star formation history in clus¬ 
ter populations. It can also be followed as a function of 
redshift (the Butcher-Oemler effect; Butcher & Oemler 
1984), a key parameter for numerous galaxy formation 
and evolution scenarios. 

The examination of the CMR in our sample begins 
with the construction of a fiducial CMR based on a sam¬ 
ple of 1104 ellipticals in the local Universe (Schombert 
& Rakos 2009b). The midline to that CMR, based on 
narrow band colors between 4100A and 5500A, is well 
described by the 12 Gyr old multi-metallicity models con¬ 
structed in Schombert & Rakos (2009a). These models 
adopt the SED spectra from Bruzual & Chariot (2003) 
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Figure 3. NEWFIRM view of sources (green boxes) within ~ 1' of 3C 068.2. The radio source is indicated by a red arrow and red square. 
The blue ellipses are X-ray point sources found with Chandra (Wilkes et al. 2013). 


for single metallicity and age populations and convolve 
each age/metallicity with a star formation history (ini¬ 
tial burst) and chemical enrichment scenario as outlined 
in Schombert & Rakos (2009a). The result is a per¬ 
fect match to the optical CMR using only old (12 Gyrs) 
populations of varying metallicity. Interpretations of 
the CMR implying the need for younger stars (frosting 
models, see Trager et al. 2000) appear to be the con¬ 
sequence of neglecting the old, metal-poor (i.e., blue) 
component of the metallicity distribution (Schombert & 
Rakos 2009a). This has been independently confirmed 
by Bregman et al. (2006) analyzing the infrared SEDs of 
early type galaxies. 

The CMR defined by optical colors is not directly 
transferable to this project, as the JK observations for 
this sample do not use redshifted filter sets as employed 
by the Rakos & Schombert program. Thus, we must 
convolve the JK bandpasses to the rest frame model 
spectra. The three cluster redshifts (1.10, 1.57, 1.83) 
place the center of the J filter at rest frames 6,000A, 
4,903A, 4,452A and the K filter at 10,57lA, 8,775A, 
7,845A, respectively. This places the observed colors be¬ 
tween the Johnson B and / range, redward of the 4000A 
break. As the typical elliptical has a spectral intensity 
peak at 4,600A (before a sharp drop off at the 4,000A 
break and the metal line rich region between 4,200A and 
4,800A), our expectation is that the observed J — K col¬ 
ors will increase (redden) from z =1.10 to 1.57, then 
shift blueward at 2 =1.83 as the J filter passes over the 
4,600A bump, as seen in Fig 7 for both a young and 


old galaxy. Indeed, that is what is seen in the mean 
cluster colors of J - K s = 1.47 ± 0.04 (CL J0442+0202), 
J — K s = 1.59 + 0.06 (3C 068.2) and J-K s = 1.46 + 0.06 
(MS 1426.9+1052), which is the uncertainty-weighted 
mean of all candidate galaxies. 

Comparison to the 12 Gyr old model CMR of 
Schombert & Rakos (2009b) is made by redshifting the 
model spectra to the redshift of the sample clusters. The 
result is shown in Fig 8 (red lines). These CMRs typically 
lie redward of the observed colors. This is, of course, due 
to the fact that the model CMR is produced for 12 Gyr 
populations. The age of the universe at the redshifts of 
the clusters will be ~5.5, 4 and 3.5 Gyrs, respectively. 

Altering the CMR for age simply requires recalculating 
the population model by halting at a younger age. The 
only color changes are due to age effects (shifts in the 
turnoff point to the blue) as metallicity effects are minor 
since single burst models have a rapid enrichment sce¬ 
nario (e.g., the typical early-type galaxy reaches 90% fi¬ 
nal metallicity in 3 Gyrs due to the rapid onset of galactic 
winds). Recalculating the models (and chemical enrich¬ 
ment) to those cluster ages produces the blue CMR lines 
in Fig 8. While the correction due to age is relatively 
small, the age corrected CMR’s are in better agreement 
with the cluster data. The close match between the red 
sequence and the model CMR’s suggests that these rich 
cluster ellipticals are effectively coeval with an early for¬ 
mation epoch. 

The formation redshifts of these cluster candidates are 
difficult to constrain within the resolution of our data. 
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Figure 4. A color-magnitude diagram of galaxy sources surrounding 3C 068.2, which is represented by a red square. The radio source 
is not the brightest cluster galaxy, being located within the blue cloud with a measured star-formation rate of ~ 390 Mq. yr 1 (Barthel 
et al. 2012). The blue, downward triangles are Chandra point sources (Wilkes et al. 2013), and are the two brightest galaxies on the flat 
red sequence (J — K s ~ 2). The larger the size of the symbol, the more significant overdensity the galaxy resides in. The largest symbol 
corresponds to a galaxy within the area corresponding to at least a lOcr overdensity. Within this region, ~4 galaxies are expected to be 
unassociated with the cluster. 


For instance, using Zf = 5 within our simulations pro¬ 
duces red sequences that are offset from the blue lines 
in Fig 8 by A(J - K s ) = -0.08,-0.09,-0.12 for red- 
shifts oi z = 1.1,1.57,1.83, respectively. As discussed 
in Section 4, the K' photometry of 10 galaxies in CL 
J0442+0202 from Stern et al. (2003) show a bias when 
compared to our K s photometry of K' — K s = 0.32±0.16. 
For these ten galaxies, this bias could imply a bluer stel¬ 
lar population, and thus a smaller Zf. Within the un¬ 
certainties of our data, we cannot formally fit a Zf to 
the cluster candidates, but the early-type galaxies are 
consistent with Zf > 5. 

Constraining the formation age of the galaxy in high 
redshift is challenging because of the inverse relation be¬ 
tween redshift and cosmic age. For example, Eisenhardt 
et al. (2008) found that for a sample of IRAC-selected 
clusters at redshifts 2 > 1, / — [3.6] colors favored Zf > 3, 
but could also be consistent with Zf = 30 (see their Fig 
19) with significant scatter in clusters at redshifts above 
z > 1.5. Our results are consistent, with at least a few 
galaxies appearing on the red sequence as old populations 
for their redshift. 

Absolute M 55 oo AB magnitudes were computed by 
convolving the redshifted J filter that each cluster is 
observed in with V. As illustrated in Fig 7, the fil¬ 
ters are compressed at higher redshift, which lowers the 
flux through the filter, so a Af-correction of -0.23, -0.50, 
and -0.75 magnitudes was applied to CL J0442+0202, 3C 
068.2, and MS 1426.9+1052, respectively. 


Clear, early evidence of evolution in the stellar pop¬ 
ulations of cluster galaxies is the Butcher-Oemler effect 
(Butcher & Oemler 1984). Simply stated, the Butcher- 
Oemler effect is a sharp increase in the fraction of blue 
galaxies in a cluster as a function of redshift. Rakos & 
Schombert (1995) studied 509 galaxies in 17 clusters be¬ 
tween redshifts of 0 and 1, and found a linearly increasing 
/b from present-day value of 5% to 80% at a redshift of 1. 
That study utilized Stromgren photometry to determine 
a color separation between nearby star-forming galaxies 
and ellipticals (anything blueward of bz — yz = 0.20 or 
vz — yz = 0.40), and then used customized filter sets to 
match those rest frame colors at the respective redshifts 
of each cluster. Any source more blue than this color 
separation is considered a member of the blue fraction 
for the cluster. 

To calculate the blue fraction using J, K s filters, we 
computed the probability distribution for the color of 
each galaxy, and the fraction that falls 0.2 mags bluer 
than the Zf = 10 CMR models of the red sequence (blue 
lines in Fig 8) are considered to be blue galaxies. The 0.2 
magnitude offset is adopted from the example of Butcher 
& Oemler (1984), and we assumed Gaussian errors for 
the probability distribution. This results in /b values of 
40+14% in CL J0442+0202, 41+18% for 3C 068.2, and 
49+17% within MS 1426.9+1052. The /b are minimally 
reduced (~ 5%) if the blue fraction is simply defined as 
the number of galaxies that are 0.2 mags bluer than the 
Zf = 10 CMR models. Errors are based on an assumed 
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Figure 5. K s image around MS 1426.9+1052 with the SQIID camera. The large, saturated object on the left-hand side of the image is 
a K s = 9.17 mag star in the 2MASS Point-Source Catalog. The radio source is indicated by a red arrow and box, while the green ellipse 
indicates a galaxy with a photometric redshift of z = 1.265 (Richards et al. 2009). The green boxes are the remainder of galaxy sources 
found within the field. The image size is 120" x 120". 


Poisson distribution (y/N a bj ), as in Rakos & Schombert 
(1995). A summary of the fg values can be found in 
Table 6. 

The fg values are less than expected from the rapid 
increase in fg from Rakos & Schombert (1995) (see their 
Figure 4). However, this can be interpreted as a conse¬ 
quence of the cosmic star formation rate as a function 
of redshift, in which galaxies at z > 1 have nearly con¬ 
stant SFR (Madau et al. 1996), while the rate steadily 
declines from z = 1 — > 0. Updated versions of the 
“Madau diagram” confirm this trend out to z ~ 3 (Perez- 
Gonzalez et al. 2005). Our values of the blue fraction 
also agree with the more intensive spectroscopic GMASS 
study (Cassata et al. 2008) which finds fg of 50% at red- 
shifts of 2. Indeed, Newman et al. (2014) find a cluster 
at z = 1.8 that hosts 15/19 spectroscopically confirmed, 
quiescent members, further corroborating blue fractions 
that are sometimes small at this epoch. These differences 
could point to large variations in stages of cluster assem¬ 
bly and evolution in the Universe at these high redshifts. 

If the sample is restricted to the regions of maximum 
overdensity within each cluster (e.g. > 10cr, discussed 
in Section 3), the blue fraction decreases in the two low 
redshift targets but increases for MS 1426.9+1052. CL 
J0442+0202 has a relatively minor decrease to fg = 
38 ± 29%, while 3C 068.2 at z = 1.57 decreases from 
f B = 41 ± 18% -+ 27±1|%. MS 1426.9+1052 at z = 1.83 
increases its blue fraction within the region of largest 
overdensity by 11% to fg = 60 ± 33%. The large vari¬ 
ability between candidate clusters at a range of redshifts 


is an interesting test of the evolution for cluster galax¬ 
ies, as the innermost regions should have had the most 
dynamical interaction. 

We also define the faint red fraction of galaxies (ff r ) as 
those galaxies that are 1.5 M 5500 AB magnitudes fainter 
than the spectroscopically confirmed BCG for each can¬ 
didate cluster and 0.2 mags redder than the CMR model 
red sequence line used to compute the blue fraction, once 
again using the probability distribution with assumed 
Gaussian errors for each galaxy. CL J0442+0202 has a 
ff r = 25 ± 14% for the entire sample of galaxies mea¬ 
sured, and ff r = 27±27% within the region of > lOcr 
overdensity. At the redshift of z = 1.57, 3C 068.2 has a 
ff r = 26 ± 18% overall, but a 22 % increase within the 
central overdensity to ff r = 48 ± 33%, which corrobo¬ 
rates the results from its /& measurement. Our highest 
redshift candidate MS 1426.9+1052 has a comparatively 
low red fraction within the overdense region of a mere 
ff r = 10;/^%, with ffr = 24 ± 17% for the entire sam¬ 
ple. These faint, red galaxies are presumably reactants in 
future dry mergers that would be hierarchically accreted 
to form galaxies on the red sequence. The blue and faint 
red fractions are listed in Table 6 . 

Simple scenarios to build up the red sequence via the 
quenching and merging of blue galaxies, followed by 
dry mergers have been proposed by Faber et al. (2007). 
These would require a significant population of interact¬ 
ing blue disk galaxies to transform relatively weak CMRs 
at high redshift into the developed red sequence features 
seen in rich, nearby clusters (see Figs 3, 4 in Eisenhardt 
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Figure 6. The color-magnitude diagram for galaxies surrounding the central radio source (red square). The green upwards triangle is a 
probable interloping galaxy with a photometric redshift oi z = 1.265 (Richards et al. 2009), when compared to the redshift of the radio 
source at z = 1.83. The largest symbols represent galaxies found within the area corresponding to a 10cr overdensity, of which roughly 4 
are expected to be field galaxies. 


et al. 2007). Papovich et al. (2012) conclude that a se¬ 
ries of dry mergers is the dominant mechanism in the 
growth of quiescent galaxies based on an in-depth analy¬ 
sis of a cluster at z = 1.62. The processes that shape the 
growth of even the most massive galaxies do not seem to 
be resolved. Some recent studies conclude that ongoing 
mergers after z ~ 1.5 are necessary to match the present- 
day BCGs (Ferreras et al. 2012; Fassbender et al. 2014), 
while Fritz et al. (2014) suggest that major mergers are 
of minimal importance. 

The extent to which there is a well defined red sequence 
CMR at high redshift impacts our models of galaxy for¬ 
mation. CMRs out to redshifts of 1 are well established 
(see Mei et al. 2009 and references therein) as well as 
in some protoclusters at redshifts of up to z = 4 (Steidel 
et al. 2005; Overzier et al. 2008). The stellar populations 
of galaxies on the red sequence at z = 1.5, which match 
passively evolving color evolution models, is consistent 
with the massive, bright end of the CMR being in place 
by a few Gyrs after the dawn of time. This suggests that 
at least some massive galaxies experience the epoch of 
primary star formation at quite high redshift. Our re¬ 
sults are consistent with the findings of Eisenhardt et al. 
(2008) (see their Figs 17,19). 

This is in agreement with the observation that proto¬ 
cluster galaxies are, in general, older than field galaxies 
of the same redshift (Peter et al. 2007), as if high density 
environments accelerate galaxy color and morphological 
evolution. Larger intrinsic scatter in the CMR with red¬ 
shift has been interpreted as due to age effects (Mei et al. 


2009); however, it becomes increasingly difficult to sepa¬ 
rate passive age changes in ellipticals with increasing star 
formation in the blue cloud. The separation between the 
red sequence and the blue cloud becomes problematic 
inside of 2 Gyrs from the onset of galaxy formation. 

The red sequence is nominally identified with ellipticals 
in the cluster core region, for they have the reddest colors 
out to high redshifts (Mei et al. 2009). This is supported 
by the fact that the slope and zeropoint of the CMR 
is unchanged (aside from age corrections) between our 
sample (at a mean redshift of 1.5) and the present-day. 
The galaxies producing an increasing blue fraction with 
redshift are perhaps the cluster proto-SO population, for 
they have the largest color scatter in present-day clus¬ 
ters and would match the “quenching” scenario, where 
infalling galaxies interact with the cluster environment 
to halt star formation (Rakos & Schombert 1995; Faber 
et al. 2007). 

While the results presented herein are similar to other 
studies based on spectroscopic and space-imaging mor¬ 
phology, we achieve many of the same goals with respect 
to the star formation history of cluster galaxies, but with 
a much smaller amount of telescope time. In addition, 
the identification of the high redshift clusters themselves, 
as preliminary observations before spectroscopy on large 
aperture telescopes and limited space-imaging time, ar¬ 
gues for an extensive ground-based near-IR survey of the 
environment around radio sources. The analysis of the 
photometric CMR, combined with accurate stellar pop¬ 
ulation models, provides a different window into the star 
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Figure 7. Model SED for an old (12 Gyr) galaxy in red, with the blue spectrum illustrating a young galaxy more representative of the 
expected SED at the high redshift of the cluster candidate galaxies. Overlaying the SEDs are the de-redshifted J filters for each candidate 
cluster, with V as a reference. At increasing redshifts, the observed filters narrow and measure less flux, so /^-corrections were applied to 
the computed A /5500 A B magnitudes for each source. Flux in arbitrary units. 


formation history of cluster galaxies. 

6. SUMMARY AND FUTURE WORK 

We have presented three candidate galaxy clusters 
spanning the redshift regime of z = 1.10 — 1.83 iden¬ 
tified as near-infrared galaxy overdensities surrounding 
radio-loud sources. From deep NIR imaging, color mag¬ 
nitude diagrams were constructed and compared to the 
galaxy models of Schombert & Rakos (2009b). We find 
that the ellipticals on the red sequence within each clus¬ 
ter candidate were all formed at a relatively old age. 
The data are consistent with Zf > 5 in our models, but 
we lack the resolution to formally fit a formation red¬ 
shift. These apparently old, massive red sequence galax¬ 
ies, while clearly present, are a minority of the cluster 
population. There is a slowly increasing blue fraction of 
galaxies within each cluster candidate as a function of 
redshift, from fn = 40% at z = 1.11 up to fs = 48% at 
z = 1.83. 

The candidate clusters are prime targets for further 
study. Use of the Hubble Space Telescope’s superior 
imaging capabilities, compared to ground-based NIR im¬ 
ages, could provide galaxy morphology and mass esti¬ 
mates (Cassata et al. 2013), along with measurements 
of individual SFRs (Lee et al. 2014b), even at these red- 
shifts. The two candidates with larger redshifts, 3C 068.2 
and MS 1426.9+1052, would benefit from multi-object 
spectroscopy, primarily to establish cluster membership, 
while secondarily providing insight into the types of 
galaxies inhabiting proto-clusters (e.g. emission line sig¬ 
natures), as Stern et al. (2003) did with CL J0442+0202. 

In the X-ray, Chandra or XMM -Newton could poten¬ 


tially be used to confirm these objects as clusters via 
their ICM emission, as well as providing a mass-estimate 
for the system. Similarly, one could use CARMA to mea¬ 
sure the SZ signal, which will accomplish similar goals as 
the aforementioned X-ray observations, but with a po¬ 
tentially more modest investment of telescope time. 

Most importantly, a confirmation of the presence of a 
physical structure with spectroscopic redshifts in any of 
these potential clusters will add to the small but growing 
number of high-redsliift clusters. The number density as 
a function of redshift and an estimate of the structures 
mass can provide constraints on cosmological parameters 
and offer comparisons to simulations of large-scale struc¬ 
ture. The leverage to constrain models improves with 
increasing mass and redshift. Further investigations into 
the properties of the constituent galaxies (morphology, 
colors, spectral signatures) will be facilitated by first es¬ 
tablishing a quasi —virialized structure. 
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Figure 8. Each candidate clusters CMR, with absolute M5500 AB magnitudes computed with the assumption that all sources are at the 
same redshift as the radio-loud source. As in the case of CL J0442+0202, there will be foreground and background interlopers within these 
CMRs. A /^-correction has been applied to each cluster candidate based on the SED template age and redshift with respect to the M5500 
magnitude system (Fig 7). The red CMR line represents the red sequence from the 12 Gyr galaxy models of Schombert & Rakos (2009b), 
redshifted to the respective value. The blue lines are recomputed SED models for the assumed ages of each candidate cluster, which are 
listed in the parentheses next to the redshift of the radio source. 
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Table 4 

3C 068.2 Sources 


Table 5 

MS 1426 Sources 


Galaxy 

ARA 

(”) 

ADEC 

(”) 

Ar 

(”) 

K s 

(mag) 

J — K s 

(mag) 

i 

0.000 

0.000 

0.000 

18.33i0.17 

1.60i0.24 

2 

-13.247 

39.826 

41.971 

18.61i0.19 

1.19i0.24 

3 

4.572 

34.854 

35.153 

17.94i0.15 

1.93i0.21 

4 a 

-7.548 

27.028 

28.063 

16.30i0.10 

1.99i0.13 

5 a 

-12.031 

21.261 

24.429 

17.55i0.13 

2.09i0.19 

6 

31.568 

20.292 

37.528 

19.26i0.25 

0.92i0.32 

7 

53.111 

5.645 

53.410 

18.84i0.21 

0.94i0.26 

8 

-14.288 

0.963 

14.321 

18.68i0.19 

0.78i0.23 

9 

-46.977 

-6.956 

47.489 

18.34i0.17 

0.98i0.21 

10 

-12.658 

-13.683 

18.640 

18.84i0.21 

1.04i0.26 

11 

12.944 

-43.318 

45.211 

19.28i0.25 

1.42i0.35 

12 

-1.630 

-20.314 

20.379 

19.42i0.27 

1.89i0.43 

13 

16.984 

-29.481 

34.023 

20.17i0.40 

0.59i0.47 

14 

12.224 

4.321 

12.966 

19.98i0.38 

4.06i2.14 

15 

-0.121 

-8.203 

8.204 

19.41i0.29 

2.94i0.80 

16 

-6.447 

1.032 

6.529 

19.67i0.32 

2.28i0.61 

17 

20.055 

-23.083 

30.578 

19.67i0.32 

2.64i0.69 

18 

-22.220 

-25.015 

33.459 

20.02i0.38 

3.30il.25 

19 

-30.541 

-30.224 

42.968 

19.76i0.31 

0.52i0.37 

20 

28.521 

-41.641 

50.472 

19.99i0.37 

1.62i0.55 

21 

35.474 

-39.058 

52.763 

20.89i0.64 

2.10il.21 

22 

52.244 

-32.845 

61.710 

19.99i0.36 

0.99i0.46 

23 

19.391 

4.007 

19.801 

19.91i0.35 

1.50i0.50 

24 

14.409 

-18.372 

23.348 

19.23i0.26 

1.63i0.36 

25 

-13.351 

-20.148 

24.170 

19.38i0.27 

1.96i0.43 

26 

28.742 

-4.167 

29.043 

19.28i0.26 

3.67i0.90 

27 

22.496 

24.888 

33.549 

19.28i0.26 

2.36i0.47 

28 

0.204 

22.287 

22.288 

18.95i0.22 

2.48i0.41 

29 

25.898 

-2.416 

26.010 

19.37i0.27 

2.15i0.45 

30 

19.655 

23.165 

30.380 

19.58i0.30 

2.24i0.53 

31 

26.082 

-34.725 

43.430 

18.74i0.20 

1.93i0.31 


Note. — Galaxy candidates within ~ 1', which corresponds 
to ~ 0.5 Mpc at this redshift (Wright 2006), around 3C 068.2 
(Galaxy 1) located at RA=2:34:23.8, DEC=+31:34:17. 
a X-ray point source identified with Chandra (Wilkes et al. 2013). 
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Table 6 

Blue and Faint Red Fractions of Candidate Cluster Galaxies 
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Redshift 

fb 

ffr 

CL J0442+0202 

2 = 1.10 

40+14% (38+29%) 

25+14% (27^27%) 

3C 068.2 

z = 1.57 

41+18% (27^2?%) 

26+18% (48+33%) 

MS 1426.9+1052 

z = 1.83 

49+17% (60+33%) 

24+17% (10±fg%) 


Note. — The blue (/&) and faint red fractions (// r ) of galaxies for each 
cluster candidate. Values in parentheses are the fractions contained within 
the regions of greatest galaxy overdensity (> 10cr). The blue fraction is 
calculated by using our Zf = 10 CMR models for each cluster candidate 
(blue lines in Fig 8) which represents the expected red sequence at that 
redshift, and specifying the color cut as 0.2 magnitudes bluer than the 
expected red sequence, following the example of Butcher Sz Oemler (1984). 
The probability distribution of each galaxy that falls below that color cut 
is deemed part of the blue fraction. The faint red fraction utilizes the 
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than the spectroscopic BCG in each cluster. Poissonian errors ( yjN 0 bj ) 
are included. 
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